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Abstract 

In this work we give detail consideration of the possible scenario of evo- 
lution of isolated neutron stars (INSs) and determine some characteristics 
of X-ray pulsars from their spin period evolution. 

The new points of our consideration are: 
-we give additional arguments for the short time scale of the Ejector stage 
( « 10 7 - 10 s yrs ). 

-we proposed specific SPINDOWN THEOREM and give some arguments 
for its validity. Discovery of accreting INSs will means that the SPIN- 
DOWN THEOREM is true. 

-we consider firstly evolution of spin period of a NS on the Accretor stage 
and predict that its period > 5 • 10 2 sec and INSs can be observed as 
pulsating X-ray sources. 

-we modeled accretion onto an INS from the interstellar medium in the 
case of spherical symmetry for different values of the magnetic field strength, 
ambient gas density and NS's mass. The periodic sources with P from 
several minutes to several months can appear. 

-we consider new idea of stochastic acceleration of very old NSs due to 
accretion of turbulizated ISM. 

-last point means that INSs can be spin up and spin down with equal 
probability. 

-using the observed period changes for four systems: Vela X-l, GX 301- 
2, Her X-l and Cen X-3 we determined D, the 'diffusion coefficient ',- 
parameter from the Fokker-Planck equation. 

-using strong dependence of D on the velocity for Vela X-l and GX 301- 
2, systems accreting from a stellar wind, we determined the stellar wind 
velocity. For different assumptions for a turbulent velocity we obtained 
V = (660 - 1440) km ■ s _1 . 

-we also determined the specific characteristic time scales for the 'diffusion 
processes' in X-ray pulsars. It is of order of 200 sec for wind-fed pulsars 
and 1000-10000 sec for the disk accreting systems. 
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1 Introduction 



Isolated neutron stars (INSs) which are not observable as radiopulsars became 
recently to attract much interest. Treves & Colpi (1991) suggested that INSs 
accreting from interstellar medium (ISM) can be observed in UV and X-rays 
by ROSAT and several of them were observed. 

For main-sequence stars there are two the most important parameters which 
are needed for description of their evolution: the mass and the age. For NSs 
there are also two important parameters: the period and the gravimagnetic 
parameter,?/ . 

The gravimagnetic parameter was firstly introduced by Davis & Pringlc 
(1981) as: 

M m 

y = V- (1) 

There are four possible states of a NS in low density plasma: E (Ejector), P 
(Propeller), A (Accretor) and G (georotator) , depending on relations between 
four characteristic radii: Ri, radius of the light cylinder, R st , stop radius , 
Rg — ( 2 ^ M )~radius of gravitational capture and R c — (^^) ^-radius of 
corotation. 

As the result we have two critical periods: Pe and Pa, separating different 
stages of NS. If P < Pe we have Ejector , if Pe < P < Pa we have the 
Propeller stage and if P > Pa and R s t < Rg -the Accretor stage. In some 
cases it may be P > Pa but R st > Rg and accretion is not possible because 
Geo-like magnitosphere is formed. 

Usially a track of a NS can be represented as a monotonic spin-down, con- 
sequantly passing through three stages: 

E — > P — > A 
E — > P — ► G 
The stage of Georotator can take place if: 

y<VG = j2GW (2) 

On the Fig.l we present examples of the NS evolutional tracks. 

In this poster we show some results, concerning periods of NSs, their evo- 
lution, determination of the stellar wind velocities in X-ray pulsars and INSs. 
Mostly we are based on three our articles (Popov 1994, Lipunov and Popov 
1995a,b) 
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2 The Spindown theorem. 

The equation of evolution of INSs can be represented in the form of the spin- 
down equation (Lipunov 1982): 

dho_ _ _hj^_ 
dt ~ Rf [6) 

where fc t -dimensionless factor, i? t -the scale of the interaction between the mag- 
netic field of the NS and surrounding plasma. 

For the next consideration we proposed some kind of astrophysical theorem 
from the general point of view. 

THE SPINDOWN THEOREM. 

The duration of the Propeller stage under all other constant conditions is 
smaller than duration of the Ejector stage: 

t P < t E , (SDT) 

Arguments. 

I) . For the most investigated propeller regimes the spin-down coefficient, k u 
slowly depends on the spin frequency, to, and as result the duration time of the 
Propeller stage is determined by the initial (not final, as for Ejector) period of 
theNS! 

II) . The interaction between the magnetic field of the NS and accreting 
plasma in each case is more effective than the interaction between the magnetic 
field and vacuum (magnito-dipole spin-down). Realy, we have the relation: 



(Ht)e _ 2 (Ra 
(*£) P 3h \R t 



Because (Ra/Ri) 3 is a very small parameter we can expect that [3 < 1, for 
a very large range of kt ■ 

Conclusion. 

If l).and 2). are true then the duration of the Propeller stage is less or equal 
than the duration of the Ejector stage. 

If we believe in the SPINDOWN THEOREM we conclude that after lO 7 " 8 
yrs the INS in real interstellar medium comes to the Accretor regime. 

We have strong observational argument for the SPINDOWN THEOREM: 
the X-ray pulsar X Per is a very slowly rotating NS (period p=835 s) associated 
with young massive companion (age of companion ss 10 7 yrs) and has a very 
low accretion rate (about ~ 10 13-14 g / s approching to interstellar medium) . We 
can not understand the existance of this X-ray pulsar if the propeller mechanism 
is not effective! 
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3 Stochastic acceleration of INSs. 



Now we will discuss the situation with very old INSs and their periods. Let us 
consider the fluctuations of the period of the NS. Suppose that the change in 
the angular velocity of the accreting star occurs under the action of a random 
torque: 

here $ is the fluctuating torque due to accretion from turbulizated ISM. We can 
introduce the scalar potential V (see Lipunov 1992): 

% c = v^M (6) 

The distribution function among frequency, f(u>), can be recieved as a result 
of stationary solution corresponding the Fokker-Planck equation: 

f(w) = N-e X p(-V(u)/D) (7) 
here D-the diffusion coefficient. For our case we have: 

V(u) = kt ^ I u I 3 (8) 
1 ' 3GMI 11 y ' 

and 

(., 

where v t and R t -velocity and length of the turbulence. We must put R t sa R G , 
and after some calculations we obtain: 

D « 7.6 • 10- V-^X^ 2 (^) 7 (10) 

Assuming V(ui) « D we obtain the estimate of the possible period of INS 
due to stochastic acceleration: 



1/3 

/ hi i I \ 

2tt 



( v 

I 3GMID J 



In principle this estimates show that real period of INS can be about several 
hours or days. But due to the very strong dependence on the velocity,^, and 
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the turbulent velocity, vt, at the scale about Rg we cannot give more precise 
results. We only illustrate that the old INSs can be accelarated by stochastic 
angular momentum from the turbulisated ISM. In this case we can see with 
equal probability spin-up and spin-down INSs! 

Spin-up and spin-down time can be estimated as: 

Ico Ivt, 

tsu ^ t s d 



Mv t R G (2GM) 3 pooWtP 

«20^/««gp=L (^"XlO^)" 1 (12) 

Rcaly vt is a small fraction of the sound velocity, because in opposite case the 
energy of turbulence will dissipate in the form of shock waves. From observations 
it is well known that turbulence has a power spectrum: 



v t ~ l a , (13) 

(here I is the scale of the turbulent torque with turbulent velocity v t ), but exact 
value of a is unknown. Very often the turbulence spectra is represented in the 
form of the Kolmogorov spectrum: 

vt = d, (14) 

where e is a constant which characterises energy transfer. If we want to 
calculate Vt on the scale of Rg we must know vt on the characteristic scale. 
Turbulence is observed on scales 10 8 < I < 3 • 10 20 cm and for I = (50 — 150) pc 
Vt is of order of 10 6 cm/s (Ruzmaikin et al 1988). we also use another value: 
Vt = 10 5 ( 10 i' cm ) 1/3 (Kaplan & Pikelner 1979). so for R G = 4 • 10 14 cm we have: 
v t « 10 4 — 10 5 cm/s. In this case t su and t s d are of order of 200-2000 yrs. In 
dense clouds (p_ 24 w 100) v t « 3 • 10 3 - 10 4 cm/s (Canuto & Battaglia 1985) 
and t su w t s d ~ 20 — 60 yrs and NS can be observed in principle. 



3.1 Accretion onto OINSs and periodic sources 

We modeled accretion onto an INS from the interstellar medium in the case of 
spherical symmetry for different values of the magnetic field strength, ambient 
gas density and NS's mass. We tried to verify the idea that if the radius of 
corotation, R co , is less than the Alfven radius, Ra , the shell will form around 
the INS and Ra will decrease to R co and the periodic X-ray source will appear 
(see Colpi et al., 1993). 

Dependence of Ra from t in our model coincides well with the analytic for- 
mula from Colpi ct al.(1993). The periodic sources with P from several minutes 
to several months can appear (see Popov 1994). 
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4 The 'Diffusion coefficients' and the stellar wind 
velocities for X— ray binaries. 

Here we shall try to estimate the diffusion coefficients for selected systems (see 
section " Stochastic acceleration of INSs" ) and stellar wind velocities (see details 
in Lipunov & Popov, 1995a, b). At first we shall estimate D. All variables, 
except At , are known (in principle). Their values taken from Lipunov (1992) 
are shown in table 1. 

Characteristic time At for the wind-accreting systems can be determined 
from the equation: 

At « 1.7 • 10 4 a- 2 10 2(A+8 - 5) i3 7 ^/i3o^ec, (15) 

where a - a fraction of the specific angular momentum of the Kepler orbit at 
the magnitosphcric radius, ^4-noise level (see table 2) (de Kool & Anzer 1993). 
So we can write equation for D in the form: 

D = 5.55 • 10- ^o-'V 1 ^ 1 (y^y *" 3 • (16) 

Values of D for four systems are shown in table 3. 

For characteristic time in the frequency space we can write: 

* - /• - Dp4 c\v\ 

£ char — W c har/U — ~^%£2 { ' > 

We can give a physical interpretation for t c h ar for wind-fed pulsars as a 
characteristic time of the momentum transfer: 

These characteristic time scales are also shown in table 3. 
Now for D we can write: 

D = 4.38 • lO- 23 ^^- (l^y ' (W) 

As we see there is a strong dependence of D on v. So we can evaluate v (in 
this case it is the stellar wind velocity, v sw ): 

v sw = 1700 • DZ% r MUW /r Tl / W' 7 G?' 7 (y^) V7 km/s. (20) 
Values of Vsw &re shown in table 3. 



6 



Vela X-l 
GX 301-2 
Her X-l 

Cen X-3 



p, sec 
283 
696 
1.24 
4.84 



Porbi SGC 

7.7 • 10 5 
3.6- 10 6 
1.5 • 10 5 

1.8 • 10 5 



Table 1: 
L, erg/sec 

1.5 • 10 36 
10 37 
10 37 
5 • 10 37 



M /10 30 Gs ■ cm 3 
3 

120 
0.6 
5.7 



^su observe Y^S 

3000 
> 100 
3 • 10 5 
3400 



^su min i 

3000 
100 

8000 
600 



Table 2: 





A 


Lmaa; 


Vela X-l 


-9.1 


10 36.8 


GX 301-2 


-8.5 


10 37 









Table 3: 










At, sec 


7 


D, sec -3 






^char-! SCO 










(v t = 0.1 • a s ) 


Ot = a s ) 




Vela X-l 


1.5 • 10 4 


6.3 • 10 6 


8.7- 10- 24 


848 


1442 


200 


GX 301-2 


1.1 • 10 3 


2.9 • 10 6 


2- HT 21 


656 


1120 


245 


Her X-l 


4 • 10 3 


6.0 • 10 7 


4-10" 19 






960 


Cen X-3 


2 • 10 4 


9.5 • 10 5 


4.1 • 10- 17 






8500 
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More information on the topic can be found at: http:/ /xray.msu. subpolar/ 



Your comments you can send to:polar@xray.sai. msu.su or ps@sai.msu.su. 
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